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Abstract 

A quasi-two-dimensional (Q2D) electron liquid (EL) is formed at the interface 
of a semiconductor heterojunction. For an accurate characterization of the Q2D 
EL, many-body effects need to be taken into account beyond the random phase 
approximation. In this theoretical work, the self-consistent static local-field cor- 
rection known as STLS is applied for the analysis of the Q2D EL. The penetration 
of the charge distribution to the barrier-acting material is taken into consider- 
ation through a variational approach. The Coulomb from factor that describes 
the effective 2D interaction is rigorously treated. The longitudinal dielectric func- 
tion and the plasmon dispersion of the Q2D EL are presented for a wide range 
of electron and ionized acceptor densities choosing GaAs/AlGaAs as the physical 
system. Analytical expressions fitted to our results are also supplied to enable a 
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widespread use of these results. 
PACS: 71.45.Gm, 73.20.Mf 

1 INTRODUCTION 

The name electron liquid or electron gas refers to a model system formed by interacting 
dynamical electrons within a medium containing a uniformly distributed positive charge 
having no motion. The overall system is electrically neutral. As the positive background 
is rigid, it does not respond to any kind of excitation, hence, it cannot polarize, however, 
the electrons can. The three-dimensional (3D) electron liquid (EL) has been studied 
as a model system for metals |1| and the 3D positive ion liquid was proposed as a 
model astrophysical system 0. In the case of two-dimensions, the study of the two- 
dimensional (2D) EL has been driven mainly by technological advances such as silicon 
inversion layers [Q, modulation doped field effect transistors Q, intercalated graphite 
layers ||, and fractional quantum Hall effect in 2D electron systems p|. In addition to 
its technological importance the 2D EL contains rich physics due to enhanced particle 
correlations and geometrical parameters that characterize the actual realization of the 
2D system. 

The EL stayed as a problem of interest in the past few decades and intense research 
efforts lead to several advances in the field. For the 2D case the first major contribution 
was due to Stern who calculated the density-density response function of the noninteract- 
ing EL |0] , which is known in the 3D the Lindhard function. The Stern function 
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(i.e., 2D Lindhard function) immediately made the random phase approximation (RPA) 
available to 2D EL. RPA was at that time one of the most successful many-body ap- 
proaches for the EL. In 1976, Jonson || showed that for 2D EL, a many-body approach 
proposed by Singwi and coworkers || (referred to as STLS) performed remarkably better 
than RPA. We have very recently compared the 2D-STLS technique with the quantum 
Monte Carlo data of Tanatar and Ceperley |T(| and proposed analytical forms for the 



dielectric function of the ideal 2D EL based on the STLS technique ||TT . 

The knowledge of the dielectric function and the local-field correction paves the way 
for a variety of many-body related terms such as self-energy, carrier lifetime, mobility 



12|| . Connections to density-functional theory can also be established [13|. Further- 
more, the dielectric screening plays a substantial role in the characterization of other 
excitations, such as, polarons jTJJ], |T5|] , |TjJ. 



In this work, our aim is to present an accurate and systematic characterization of 
the dielectric properties of the quasi-2D (Q2D) EL in real heterojunctions where the 
electron distribution can penetrate to both sides of the interface. The charge distribution 



is based on a variational approach proposed by Bastard ||17|| . The effective 2D electron 
interaction for this system is characterized by the Coulomb form factor. This quantity 
is treated rigorously. The dielectric function for the GaAs/AlGaAs hetero junction is 
given for wide ranges of electron and ionized acceptor densities. Throughout this work 
the term dielectric function, refers to the longitudinal dielectric function. We also fitted 
analytical expressions to our data for the efficient use of these results by other researchers. 
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To simplify the computational labour we stayed in the zero-temperature formulation 
and the so-called electrical quantum limit, where only the lowest subband along the 



confinement direction is populated (we refer to a very recent work JT8|, discussing the 
effects of higher subbands on the dielectric function). 

The paper is planned as follows: Sec. II discusses briefly the variational computation 
of the Q2D electron distribution. The effective 2D interaction of these Q2D electrons 
is treated in Sec. Ill and the modifications to the STLS technique in the Q2D case is 
contained in Sec. IV. In Sees. V and VI the dielectric function and the plasmon dispersion 
are considered respectively; all the results are given referring to GaAs/AlGaAs as the 
physical system, however, the approach is developed for a general heterojunction. In 
Sec. VII, the fitted analytical expressions for the results are presented. Following the 
conclusion section, appendices include some details on variational formulation and the 
Coulomb form factor for a Q2D system. 



2 VARIATIONAL CHARGE DISTRIBUTION FOR 
A HETEROJUNCTION 

The electrons from ionized donors in the barrier side of a modulation doped heterojunc- 
tion are trapped in a wedge-like well formed by a step barrier due to conduction band 
edge discontinuity on one side, and the potential due to presence of the transferred elec- 
trons and ionized acceptors on the other ||. The one-dimensional quantum confinement 
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gives the Q2D nature to the system and behaves remarkably different than ideal 2D 
and 3D systems. In handling the many-body effects in heteroj unctions we avoid some 
critical simplifications that have been used in the past such as infinite barrier height 
I4jl , [ D| , ||l6f (which is a reasonable approximation only for Si inversion layers) and no 



ionized acceptors within the channel JT9[ (which is in fact not the case in practise). For 
an accurate account of the electronic distribution in heterojunctions, we use Bastard's 
variational approach that was tested previously in determining the subband energies 



The electronic wave function ^(z), within the effective mass approximation satisfies 
the one-dimensional Schrodinger equation along the confinement direction (chosen to be 
the ^-direction) 
-h 2 d 1 d 



+ U e - e (z) + U A (z) + U barrier (z) 



<*(z) = Ei*(z). (1) 



2 dz m(z) dz 

U e _ e (z) is the potential (energy) formed by the presence of the electrons, U A (z) is the 
potential due to ionized acceptors and Ub arrier (z) is a step-barrier potential: Ub arrier (z) = 
UbQ(-z), resulting from the conduction band edge mismatch of the neighbouring ma- 
terials. m(z) is the effective mass of the conduction band electrons being equal to 
in the barrier-acting material and m A in the well-acting material. Bastard proposed 
the following variational form for the lowest subband Si(z) allowing penetration to the 
barrier region, [z < 0) jlTJ] 

Me KbZ / 2 , for z < 

QU)H ~ • ( 2 ) 

N(z + z )e~ bz / 2 , for z > 



Invoking the continuity of Si(z) and m (z)4-si(z) at z = and the normalization of 
Si(z), dz\qi(z)\ 2 = 1, yields the following three equations, 

M = Nz , (3) 



TUB 



^2[l + te + ^(l + A 



b 3 

iPzZ f. ~- ( 5 ) 



N 



Bastard also set Kb = 2y2rriBUi } /h 2 and used b as the only variational parameter. 
We have observed that such a choice of is highly satisfactory in the electrical quantum 
limit. Note that M,N and z also depend on b through Eqs.@-(0). b is determined 
by minimizing the total system energy (see Appendix A for the expressions). A closed 
form representation of b is not possible, unlike the Si-inversion layer H, however, the 
minimization can easily be achieved numerically. We work in the regime where only 
the lowest subband is populated, this puts an upper limit to the 2D electron density 
above which the Fermi level crosses the first-excited subband energy. For GaAs/AlGaAs 
heterojunction our analysis is valid for the 2D electronic densities, N2D < 7x 10 11 cm -2 . 
Bastard's work JT71 can be consulted for further details. 
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3 COULOMB FORM FACTOR FOR A PENETRA- 
BLE HETERO JUNCTION 

In the 2D EL the interaction potential in reciprocal space is taken to be 2ne/q, where 
q is the wave number. This potential is obtained by taking the 2D Fourier transform 



of the 3D Coulomb interaction which is 1/R, R denoting distance in real space [ 20| . In 
fact a strictly 2D solution of Poisson's equation is proportional to — \\i{R) rather 
than 1/R and its 2D Fourier transform is proportional to 1/q 2 as in 3D EL. However, 



the — \n(R) interaction is seldom used [22] due to indication by real physical 2D systems 



that 1/R type of interaction is relevant [Z5|,||24||. For the case of a Q2D system the 



charge distribution along the third dimension modifies the effective 2D interaction from 
2 r ne/q to F(q) 2ire/q. F(q) is the Coulomb form factor describing the effect of the finite 
spread of the charge distribution along the confinement direction over a region where the 
background dielectric constant is discontinous due to different materials on both sides. 

Following the approach in the previous section we use the variational charge distri- 
bution that can leak into the barrier region and calculate the function F(q) accordingly. 
The details on F(q) are given in the Appendix B, here we state the final result, 



where, 
h = 2N 4 



1 

2b 



7 A 2? 3 2z 2 

z _|_ zz _|_ zz 



(b + q) (b + q) 2 (b + qf 



(26) s 



Azl + 6z 2 + Az 



(b + q) (b + q) 2 (b + q) 3 



(2b) 3 \ (b + q) (b + q) 2 (b + q) 3 \ ' (2b)* \(b + q) ' (b + q) 2 \ ' (26) 5 (6 + g) 



I 2 = AT 



+ 



2z 



+ 



(6 + g) (6 + g) 2 (b + q) 3 
M 4 



M 4 



(«& + <?) 2 

M 2 N 2 
(«6 + q) 



2z 



(b + q) (b + q) 2 (b + qf 



(8) 

(9) 
(10) 

(11) 



In Eq. (|) 6,4 and €3 are the background dielectric constants of the well-acting and 
the barrier- acting materials respectively. The bare electron-electron interation potential 
energy for this Q2D system becomes 

(12) 



uQ 2B {q) = ^ F{g) 



eq 

where e = (e^ + €b)/2 and q is the 2D wave number associated with the spatial variation 
along the 2D sheet. 

The terms containing I 2 and I± in Eq. (|6|) represent the image interaction resulting 
from the different permittivities on both sides. Their effects decrease when the permit- 
tivity contrast diminishes; an example is the GaAs/AlGaAs system considered in Fig. 1 
for two different electron densities (see the following section for the material parameters 
used). The Coulomb form factor becomes more important in high electron densities 
(see Fig. 1) where the in-plane particle separation is comparable to the extension of the 
charge distribution along the confinement direction. The expression for F(q) in Eq. (||) 
will especially be useful for heteroj unctions with a high permittivity difference and a 



low barrier height. 



4 Q2D STLS 

The STLS technique in 2D has been discussed in the literature, and we refer, for instance, 
to Jonson's pioneering paper ||. In going from 2D to Q2D the only modification (within 
the electrical quantum limit) is the replacement of the 2D Coulomb interaction energy 
by the effective 2D interaction due to finite extension of the charge distribution along the 
confinement direction. The exchange and correlation hole associated with each electron 
in the system is described by the local-field correction, G(q). This function in the case 
of Q2D STLS reads 

G ^-II^W)l£ ll - s(W - m ' (13) 

where the subscript n is used in this equation and in the rest of the text to denote wave 
numbers normalized to the Fermi wave number kp (i.e., q n = q/kp etc.). In Eq. ([13]) 
F is the Coulomb form factor and S is the static structure factor. The latter con- 
tains contributions from plasmons and electron-hole pairs and is related to the dielectric 
function through the fluctuation-dissipation theorem. The dielectric function, in turn, 
depends on the local-field correction (see Sec. V). The computational task involves the 
self-consistent solution of these three coupled nonlinear integral equations. A change 
of variables leads to a substantial improvement in the execution speed of the STLS 
algorithm. Using t n = p n — q n in Eq. ( |13"D leads to 



G Q2D (q) = r dt n t n [1-S(t)) f #F(Vl + a2 + 2acos0) 7=^=^ -(14) 

7iF(q)Jo Jo \ v / Vl + a 2 + 2acos0 

where a = t n /q n . 



In Fig. 2 we present the self-consistent STLS G® 2D (q) results for a wide range of 
electronic densities given in terms of r s . r s is the effective interparticle spacing defined 



as r s = 1/ a* B \JnN2D where N 2 d is the 2D electronic density and a* B is the effective Bohr 
radius given by a* B = ^77^2, e is the background average static dielectric constant 
and m*m is the effective mass of the electrons considered, with mo being the free 
electron mass. We consider GaAs/AlGaAs heterojunction as the physical system with 
the parameters = 0.07mo, m B = 0.088m , €a = 13, e# = 12.1 and £/& = 0.3 eV 
(corresponding to an Al mole fraction of 0.3) which were used by Stern and Das Sarma 



2q| . For a* B we used e = 12.55 and m* = 0.07, giving a* B = 9.49 nm. The conduction 
band offset, £/& was measured by some groups to be around 0.225 eV (in contrast to 
0.3 eV) PB| . We have observed that our results are not sensitive to the deviation of £/& 
in this range. In Fig. 2 the interval r s = 0.8 — 20 is shown with an ionized acceptor 
density of N dep i = 0.46 x 10 11 cm~ 2 . For r s < 0.8 the higher subbands start to be 
populated which was not taken into our analysis. 



For the 2D EL, STLS G(q) becomes proportional to q as q — >• [IT] whereas in 3D 
case it is proportional to q 2 [IJ. In the Q2D case we observe that (see Fig. 2) for low 
r s values small-g behavior is close to quadratic and as r s increases this behavior goes 
towards a linear one indicating an approach to a 2D character. 

Gold and Calmels also reported their results on G® 2D (q) for GaAs/AlGaAs het- 



erostructure fll9"|. Their treatment is based on STLS but with essential discrepancies 



compared to ours. They imposed the local-field correction for 2D and Q2D to be of the 
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form 

n Q2D ( n_ 2/3 1.402 X 

GC {X) ~ Ts [2.644CS,(r.) + x*CUr s )}^ (15j 
where x = _ 1 1 / 3 and the coefficients Ci 2 and C 2 2 were tabulated |Tj|. They assumed 
no penetration to barrier region in the Coulomb form factor and also neglected the 
presence of ionized acceptors in the well-acting region. Especially, the form used in 
Eq. ( |15|) enabled them to reduce the computational effort appreciably, however, their 



results are in strong disagreement with ours for r s > 1 and q ~ 2k p both in 2D fll 
and Q2D as can be seen in Fig. 3. The form in Eq. (p~5|) cannot accommodate the full 
STLS G(q) leading to a poor dielectric function and screening properties. The ionized 
acceptors in the well region play a primary role and need to be included in the treatment. 



5 DIELECTRIC FUNCTION 

The function of practical importance is the wave number- and frequency-dependent 

(longitudinal) dielectric function, e(q,u) that not only determines the response to a 

weak external perturbation but also possesses information on the many-body dynamics 

of the system. With the knowledge of the local-field correction, e(q, u) is given as 

f Q2D , , . _ 1 - U Q2D (<1) 7T°(q,uj) [1 - gg°(g)j 
l + W 2D {q) 7r U (g,cu) GQ 2D {q) 

where 7T (q, u) is the 2D zeroth-order polarization insertion, the Stern function 0, [11 



Apart from 7T^, 2D and Q2D quantities behave differently. This is illustrated in Fig. 4 
showing inverse static dielectric function, e~ 1 (q,0) within RPA and STLS for both 2D 
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and Q2D cases. To assess the effect of penetration of the charge distribution into the 
barrier region, we compare Ub = 0.1 eV case with Ub — > oo in Fig. 5 at r s = 0.8. It 
is observed that for GaAs/AlGaAs-like heteroj unctions, this penetration has a minor 
effect on the static dielectric function. In Fig. 6 the inverse static dielectric function 
of GaAs/AlGaAs heteroj unction is plotted in the density range r s = 0.8 — 20 and 
for Nd ep i = 0.46 x 10 n cm -2 . Notably, the GaAs/AlGaAs heterostructure shows an 
overscreening effect (i.e., e < 0) for r s > 3. The onset of overscreening shifts to higher 



electron densities for the strictly 2D case [ II], due to enhanced particle correlations 



in lower dimensions. As an interesting consequence, the negative dielectric function 



suggests a negative compressibility of the Q2D EL |12j and in fact, recently this has 
been experimentally observed on a GaAs quantum well structure |f2~7f . 

We would like to include some necessary remarks about this dielectric function. The 
expression in Eq. (ffB) only gives the Q2D EL dielectric function. The total screened 
electron-electron interaction is 
Ittp 2 1 

UTM = F(q)— (17) 

The dielectric responses of the polar lattice and the valence electrons are contained in 
the average background dielectric constant e. Here we have used the static dielectric 
constant (see, for instance, our definition of a* B in Sec. IV), hence, it is assumed that the 
polar lattice can follow the external excitations. Obviously this limits the validity range 
of this work to ui <C ujto, with uoto being the transverse optical phonon frequency. This 
limitation is relaxed if the background lattice does not have a polar character. Hence, 
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for the particular system that we are considering , the dielectric function is expected to 
be valid up to about 1 THz. In principle, however, the static nature of the local field 
correction of the STLS technique can further limit this upper frequency. 

Finally, the dielectric function given by Eq. (|T6|) takes into account the polarization of 
the electrons in the lowest subband. Even though the presently available experiments on 



GaAs/AlGaAs systems mainly fall into this regime [pq] , fl29fl , the technological trend aims 
to populate the higher subbands to increase the amount of current carried in modulation 
doped field effect transistors by using different materials such as InGaAs/InAlAs 
When the higher subbands are occupied the dielectric function should necessarily be a 
tensor of the form eij(q, u), where i = j terms account for the intrasubband polarizations 
and % 7^ j terms represent intersubband couplings. To assess the performance of the 
presented approach regarding the electrical quantum limit, we extended the variational 
wave function technique to include lowest two subbands and determined the subband 
populations by invoking self-consistency between Poisson and Schrodinger equations. In 
Fig. 7 we show the charge distributions along the confinement direction for a density 
of 1 x 10 12 cm~ 2 . The solid curve represents the correct charge distribution containing 
contributions from the lowest and first-excited subbands. The dashed curve, on the 
other hand, sticks to the electrical quantum limit which actually breaks down beyond 
N2D = 7 x 10 11 cm~ 2 . It is important to observe that the difference between the two 
curves is quite marginal. This is simply because the percentage of the first-excited 
subband electrons is 4.7% at this density. 
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6 PLASMON DISPERSION 

The elementary excitations in electron liquids are electron-hole pair creations and col- 



lective excitations knowns as plasmons |[30|| . The latter can be characterized with the 



knowledge of the wave number and frequency-dependent dielectric function, e(q,u). 
Particularly, the plasmon dispersion relation, u p (q) is available through the zeros of the 
dielectric function; 

€(q,u p (q))=0. (18) 

Inserting the expression for e(q,u>) from Eq. ([16]) leads to the following closed form 
expression for the plasmon dispersion 



Vp{q) = - \ q 2 n + _ 2 , — , (19) 



z 2 + 2z' 



where 

q n 



V2r s F{q)\l-G^ D {q)} 

and 



(20) 



Vp{q) _J^__^ (21) 

which is valid in the range [0, q n ,max] where q n , ma x satisfies v v {q n ^ max ) = q n ^ max + q 2 n ma J1 
and outside this region plasmons dissociate to electron-hole pairs so that collective ex- 
citations are no longer long-lived. The Eq. (|T9| ) reduces to the ideal 2D result [31| when 



F(q) — > 1. Fig. 8 shows the plasmon dispersion for GaAs/AlGaAs heterostructure with 
Ndepi = 0.46 x 10 11 cm -2 and for several r s values. Even though the plasmon dispersion 
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can be experimentally probed, such as, through far infrared spectroscopy [Eg], the avail- 



able experimental results pertain to high electronic densities and small wave numbers 
(q < kp). Therefore the effects of the local-field correction have not yet been verified. 



In this section we present our fitted expressions to G® 2D (q) and F(q) applicable to 
GaAs/AlGaAs hetero junction in the density range r s = 0.8 — 20. As a fit to G® 2D (q) 
(shown in Fig. 2 by solid lines), we tried a simple form containing three fitting parame- 
ters, 



where A, B and C are the fitting parameters. The optimized values are tabulated in 
Table I for N^epi = 0.46 x 10 11 cm -2 . The third parameter, C is introduced based on our 
observations on the long- wavelength behavior of G® 2D (q) in Sec. IV. In ideal 2D, C was 
equal to one and in 3D case C was equal to two. Optimized C values in Table I show 
this interpolation between r s = 0.8 to 5, but then this trend is lost to enable a good fit 
for the whole q values. The fitted expressions are plotted in Fig. 2 by the dotted lines. 
To assess the quality of the fitting we use the following error estimate between a target 
vector, T(i) and the fitted vector Tf it (i): 



7 ANALYTICAL EXPRESSIONS 




(22) 



1 



N 



T(i)-T fu {i) 
T(i) 



error (%) 



N 



E 



100. 



(23) 



i=l 



Accordingly the deviation of the fitting in Fig. 2 is less than 2.5 %. 
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The Coulomb form factor, F(q) also requires a laborious work for GaAs/AlGaAs 
system. This function can be fitted by a simple expression 

F '«^ = irk- (24) 

containing a single fitting parameter D which is tabulated in Table I for the same N^epi 
value. 

The knowledge of G ( f^ > {q) and Ff it (q) is sufficient for representing the dielectric 
function (see Eq. (0)). The performance of fitting for e _1 (g, 0) is available from Fig. 6 
(shown by dotted lines) where the error, using the estimate in Eq. ( p3|) is less than 1%. 
Similarly in Fig. 8 the plasmon dispersions with the use of the fitted forms are shown in 
dashed lines, the fitting error being much less than 0.1%. 

We have observed that taking the barrier height £/& = 0.225 eV does not significantly 
affect the parameters A, B, C, and D. However, Nd ep i takes an important part in both 
G(q) and F(q), so we repeated the self-consistent Q2D STLS technique for N dep i = 
0.146, 1.47, 4.69 x 10 11 cm -2 and performed again fittings. Rather than specifying these 
results in tabular form we present below fitted functions of r s for A,B,C, and D. 

A m = 1.02 



1 - ai r? e- a3r °\ , (25) 
B fit = b x ln(6 2 r s ) + b 3 , (26) 



C fit = 0.42r; Cl + 1.03r s u ^, (27) 

^ = ^T7F' (28) 
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the constant parameters contained in these expressions are tabulated in Table II for the 
considered range of Na ep i values. With the expressions in Eqs. (pSj)-(p^), inverse static 
dielectric function can be generated to an accuracy of about 1 % except r s = 2 case 
having an error about 9 %. Similarly with these equations plasmon dispersion can be 
recovered to an error much less than 0.1 %. 



8 CONCLUSION 

The dielectric properties of the Q2D EL in a heterostructure are studied and the be- 



haviour is seen to be remarkably different than the strictly 2D EL |fTI ]. The analysis is 
rigorous with the only simplifications being the electrical quantum limit and the zero- 
temperature formalism. These simplifications can also be relaxed at the expense of 
computational complexity. The leakage of the charge distribution to the barrier region 
is included in the analysis through a variational approach. The full form of the Coulomb 
form factor applicable to a general heterostructure is presented. For the GaAs/AlGaAs 
system, the image terms have been observed to have a marginal role. A sizeable contri- 
bution will be encountered in the case of heterostructures built up of materials having 
a high dielectric constant contrast and a low conduction band offset. The dielectric 
function and the plasmon dispersion of the Q2D EL are characterized using the STLS 
many-body approach that leads to substantial improvement over the conventional RPA. 



Unfortunately, in contrast to the ideal 2D case ||10||, ||32|| , quantum Monte Carlo simula 



tions are not available, to compare our results, for the Q2D EL; the present experimental 
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data cannot cover the regime where the RPA breaks down (i.e., r s > 1 and q ~ k F ). Our 
analysis extends to a wide range of electron and ionized acceptor densities. To the best 
of our knowledge this work forms the most elaborate study of the screening properties 
of the Q2D EL. Our results are supplemented with analytical expressions fitted to our 
data. We have presented the expressions for the local-field correction, Coulomb form 
factor, the dielectric function and the plasmon dispersion of the Q2D EL. With this 
information, the self-energy, quasiparticle lifetime, mobility etc., can also be obtained; 
polarons in Q2D systems can be studied with the inclusion of electron-electron screening. 
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APPENDIX A: TOTAL SYSTEM ENERGY IN THE 



VARIATIONAL APPROACH 



In this section, for completeness we include the expression for the total system energy of 



a heteroj unction in Bastard's variational approach JL7). The ground-state expectation 
of the total system energy per electron is 



1 



(Er OT (6)) = (T(6)> + -(U e - e (b)) + (U A (b)) + (U bamer (b)), 



where (T(b)) is the kinetic energy term given by 



A/f 2 «- /V 2 



(U e - e (b)) is the average electron-electron interaction potential, 



(Ue-e(b)) 



87r nr 

—N 2D 
e 



/33 25te n 176 2 z 2 



iV 2 M 2 



^ 2 6 2 + 2te + 2) 



Ry, 



where (£/a(&)) is the average electron- ionized acceptor interaction potential, 

Ry, 



j7T 



(U A (b)) = —N depl 



+ M2 



b 4 



6 



where (Ub arrier (b)) is the average potential energy due to step-barrier, 



(Ubarrier (^)) 



K 6 



Ry, 



(29) 



(30) 



(31) 



(32) 



(33) 



where K b = 2Jm%Ub i ay. In above equations atomic units are used; all energies are 
in Rydbergs (lRy = moe 4 /2fr 2 ) and lengths in Bohr radius (a_B = h 2 /moe 2 ). In the 
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total system energy, the electron-electron interaction is weighted by 1/2 to avoid double- 
counting. The variational parameter b is determined by minimizing (Etot(o)) which is 
an easy task numerically as the cost function has a single minimum. 



20 



APPENDIX B: EFFECTIVE 2D COULOMB INTER- 



ACTION 

We first recall the electrostatic potential due to a point charge Q, at a distance d (along 
the z-axis) from the interface formed by two semi-infinite dielectric media with permit- 
tivities 6a and 6b ■ Solution of the Poisson's equation subject to continuity requirements 
at the interface, z = results in an electrostatic potential of the form f33fl , 



$(R,z) = < 



A \^/R2 + (d-z) 2 £ A+£B y jR2 + {d+z)' 2 



z>0 



(34) 



Q 



z<0 



This result will now be used in constructing the effective 2D Coulomb interaction energy 
between two (charge) distributions n(R, z) and n(R', z') 



U Q2D (R - it) 



I poo roo 

dz / dz' 



e A \J0 ./ii 

e 2 I f o ro 
H { I dz I dz' 



n(R, z) n(R! , z') ^tA — ^B n(R, z) n(R' , z') 

R - R'\ 2 + (z - Z') 2 e A + £B J\R_R'\2 + ( Z + Z i)2 



n{R,z) n(R',z') 



+ 



€b — €a n(R, z) n(R' , z') 



R - R'\ 2 + (z- Z') 2 e A + € B J\R - R'\2 + ( z + 



A2 



+ 



2e 2 



eA + e B Jo 



dz r V <R,z)n(S,*) 



R - R'\ 2 + (z - z') 2 



z 



° iftf n(R,z)n(R^) 
- x Jo - n R- R'\ 2 + (z- z') 2 



(35) 



The first two terms in Eq. (|35f) represent direct and image interaction of the charge 
distributions on the right side of the interface (z > 0). Third and fourth terms represent 
the same interactions for z < region. The last two terms which are in fact equal, 
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represent the direct interaction between charge distributions on opposite sides of the 



interface. 

The charged-particle distribution is n(R,z) = n(z) = \^i(z)\ 2 , where <,i(z) is given in 
Eq. (0). The 2D Fourier transform of U® 2D (R — R') is easily obtained using the result 



/ 



d 2 r f = — e-l°l« 

Vr 2 + a 2 q 



(36) 



as 



uQ2D{q) = 2_ne^ {h _ 



fiZf-j,) + (j, + !£^7 4 ) + -^f-/ 5 , (37) 
+ / <?e# V e A + e B J q(e A + e B ) 



with, 



h= dz dz'N\z + z ) 2 {z' + z ) 2 e- bz e- bz e-\ z - z \ q 
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o Jo 
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oo J — oo 



dz [° dz'M 4 e Kb(z+z ' ) e (z+z ' )q 



oo •/ — oo 



poo rO 

h= dz dz'M 2 N 2 e KbZ (z + z ) 2 e~ bz e 

Jo J-oo 



(38) 
(39) 
(40) 
(41) 
(42) 



These integrals are straightforward and the results are listed in Eqs. (|7|)- (1TT|) . The 
effective 2D interaction U® 2D (R — R') must reduce to ideal 2D \R — R'\ — > oo. 



lim U Q2D (R-R') = ^-^ 

\r-r'\^oo e A \R-R'\\ e A + e B 



1 



dz n(z) 



2 e 2 
+ 



e B \R-R'\ 



1 + 



£a + e B 



dz n(z) 



+ 



4e 2 



(e A + e B )\R-R'\Jo 
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dzn(z) / dz'n(z'). (43) 



Using J °° n(z) dz = 1 — Jl^ n(z') dz' in Eq. fl43| ) leads to the desired result, 



lim U Q2D (R-R') 



2 



(44) 



\R-R>\->oc 



R-R'\eA + e B 



As a consequence the Fourier transform gives 




2vre 2 



(45) 



qe 



so that liniq_^o F(q) = 1 as can be observed in Fig. 1. This also indicates that all of the 
interactions are properly accounted in Eq. (|37|) . 
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List of Tables 



1. Fitting parameters A, B, C, and D used in Eqs. ( P2"D and (|24]) as a function of r s 
for the characterization of the Q2D EL in a GaAs/AlGaAs heterostructure. The 
ionized acceptor density is Ndepi = 0.46 x 10 11 cm~ 2 . See Sec. IV for the other 
parameters used for GaAs/AlGaAs system. 

2. The constants used in Eqs. (|25[) -(|28| ) for different ionized acceptor densities, N dep i. 

The parameters characterizing the heterostructure are chosen suitable to the GaAs/AlGaAs 
system (see Sec. IV). 
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List of Figures 

1. The Coulomb form factor F and the effect of the image terms as a function of 
wave number q (in units of kp) for the electronic densities r s =0.8 and 20. The 
full lines apply to GaAs/AlGaAs heterostructure having 6^=13 and €£=12.1. The 
dashed lines refer to the same system but with = e#= 12.55 so that no image 
term appears. See Sec. IV for the definition of r s and other parameters used for 
the system. 

2. The local-field correction, G® 2D (q) of Q2D EL versus wave number q (in units of 
k F ) for r s values 0.8, 1, 2, 3, 4, 5, 10, 15, and 20. Solid lines: STLS and dashed 
lines: calculation using the fitted form for G® 2D (q); see Eq. (p2|) with the values 
from Table I. 

3. The comparison of the full STLS Q2D local- field correction (solid lines) with that 
of Gold and Calmels' (dashed lines) given by Eq. fllBD as a function of wave number 
q (in units of kp) for r s —l and 10. 

4. Comparison of ideal 2D and Q2D inverse static dielectric function, l/e(g, 0) as a 
function of wave number q (in units of kp) for r s =3. Solid lines: STLS, dashed 
lines: RPA. For Q2D EL, GaAs/AlGaAs heterostructure is used with N^pi = 
0.46 x 10 n cm- 2 . 

5. The effect of the barrier height, Ub on the inverse static dielectric function, l/e(q, 0) 
as a function of wave number q (in units of fcp) for r s =0.8. Solid line: Ub = 0.1 eV, 
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dashed line: Ub — > oo. Other parameters for the heterostructure are given in 
Sec. IV. 

6. The inverse static dielectric function of Q2D EL, l/e(q, 0) as a function of wave 
number q (in units of fcp) for r s values 0.8, 1, 2, 3, 4, 5, 10, 15, and 20. Solid lines: 
STLS and dashed lines: calculation using the fitted forms for G® 2D (q) given by 
Eq. (|22"D and F(q) given by Eq. ( P4[) with the values in Table I. 

7. The electron distribution along the confinement direction in arbitrary units. The 
total electron density is 1 x 10 12 cm~ 2 . Other parameters are given in Sec. IV. 
Solid line refers to the two subband populated calculation and the dashed line is 
based on the electrical quantum limit. 

8. The normalized plasmon energy (E p /Ep = 2u p ) as a function of wave number q 
(in units of kp) for r s values 1, 5, 10, and 20. Solid lines: STLS and dashed lines: 
calculation using the fitted forms for G® 2D (q) given by Eq. (p2|) and F(q) given 
by Eq. (^J) with the values in Table I. The dotted line marks the onset of the 
electron-hole continuum. 
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A B C D 



0.8 0.6243 0.4923 1.5462 1.2750 

1.0 0.6549 0.5005 1.5079 1.1542 

1.5 0.7250 0.5274 1.4342 0.9285 

2.0 0.7857 0.5519 1.3950 0.7690 

2.5 0.8380 0.5763 1.3644 0.6497 

3.0 0.8794 0.5999 1.3512 0.5571 

4.0 0.9405 0.6461 1.3274 0.4321 

5.0 0.9779 0.6855 1.3264 0.3494 

6.0 1.0012 0.7209 1.3356 0.2922 

8.0 1.0225 0.7792 1.3683 0.2197 

10 1.0294 0.8223 1.4097 0.1752 

12 1.0305 0.8597 1.4545 0.1454 

15 1.0295 0.9007 1.5014 0.1158 

20 1.0257 0.9555 1.5185 0.0863 
Table 1: 
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N dep i (cm -2 ) ai a 2 a 3 bi b 2 h c\ di d 2 

0.146 x 10 11 0.6384 0.2213 0.5555 0.3023 4.9907 -0.0435 0.6891 4.3194 2.4659 

0.46 x 10 11 0.6770 0.2794 0.6372 0.2575 2.6283 0.1623 0.7914 2.7325 1.3674 

1.47 x 10 11 0.6953 0.2302 0.6888 0.2253 1.4913 0.3043 1.1675 1.6566 0.6898 

4.69 x 10 11 0.6887 0.1802 0.7418 0.1978 0.6977 0.4195 1.5437 1.0004 0.3518 

Table 2: 
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INVERSE STATIC DIELECTRIC FUNCTION 
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